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We analyze in detail the expansion of a ID Bose gas after removing the axial confinement. We show 
that during its one-dimensional expansion the density of the Bose gas does not follow a self-similar 
solution, but on the contrary, it asymptotically approaches a Tonks-Girardeau profile. Our analysis 
is based on a nonlinear Schrodinger equation with variable nonlinearity whose validity is discussed 
for the expansion problem, by comparing with an exact Bose-Fermi mapping for the case of an 
initial Tonks-Girardeau gas. For this case, the gas is shown to expand self-similarly, with a different 
similarity law compared to the one-dimensional Thomas-Fermi condensate. 



During the last years, the achievement of Bose-Einstein 
condensation (BEC) [Q has generated an extraordinary 
interest in the physics of ultracold atomic gases. Among 
the topics related with the physics of ultracold gases, 
the issue of low-dimensionality is attracting a growing 
interest. The recent development of trapping and cool- 
ing techniques has enabled experimental realizations of 
low-dimensional gases both in one (2]-||] and two [||J|-0] 
dimensions. For the case of ultracold low-dimensional 
dilute atomic gases, it has been theoretically predicted 
that ID ||, 2D ||[l(| and even very elongated but still 
dynamically 3D gases JlT| , should present an equilibrium 
BEC with a spatially fluctuating phase. Such quasicon- 
densates have been recently observed by means of time 
of flight measurements B |. 

The impenetrable ID gas of bosons, the so-called 
Tonks-Girardeau (TG) gas, has recently deserved a spe- 
cial interest [p 13-13]. In order to accomplish this partic- 
ular regime, rather strict conditions for the temperature, 
gas density, interaction potential, and trapping potential 
must be fulfilled |^,^6| . These conditions can be achieved 
with currently available experimental techniques. Par- 
ticularly important in this sense is the recent progress in 
loading ID Bose gases in optical lattices JijJ] where the 
transversal confinement can reach 100kHz and the devel- 
opment of the Feshbach resonance techniques to modify 
the value of the s-wave scattering length |l8). Therefore, 
it is important to characterize the properties of the TG 
gas, and especially the intermediate regime between the 
quasi-lD BEC and the TG gas. 

For the TG gas with a zero-range infinitely repulsive 
interatomic potential, the bosons acquire effectively a 
fcrmionic character and the mapping between bosonic 
and fermionic wavefunctions is exact, both for homoge- 
neous []l9| and trapped gases [Q. Interestingly, the ho- 
mogeneous delta-interacting ID bosonic gas under peri- 
odic boundary conditions is analytically solvable for any 
strength of the interactions, as shown by Lieb and Lin- 



iger (LL) [gfj. There is unfortunately, to the best of 
our knowledge, no exact solution for arbitrary interac- 
tion strength in the case of trapped gases. An interest- 
ing approach was introduced in Ref. 0], where a hy- 
drodynamic formalism was shown to reproduce the sta- 
tionary properties of the TG gas. The approach of Ref. 
PH should, however, be employed carefully since it sig- 
nificantly overestimates the coherence of the system [^4| . 
Recently, the approach of Ref. |2l| was extended to the 
case of finite interactions, by employing the LL model and 
local density approximation JIB] . In Ref. |16| the density 
profile of the trapped gas was analyzed for regimes rang- 
ing from Thomas-Fermi (TF) profiles to TG. A differ- 
ent approach to the issue of finite interactions has been 



discussed in Ref. |23|, where the intermediate regime is 
considered as a mixture of a BEC and a fermionized TG 
gas. 

This Letter is devoted to the analysis of the ID expan- 
sion of Bose gas. We employ the procedure of Ref. |ll| 
to show that contrary to the case of a ID Thomas-Fermi 
condensate, the expansion is not self-similar. In fact, 
any initial density profile will asymptotically evolve dur- 
ing its expansion towards a TG-like shape. Consequently 
the ID expansion allows us to easily explore intermedi- 
ate situations between the TF and TG regimes. Addi- 
tionally, we discuss with the help of the Bose-Fermi (BF) 
mapping jl4|,[l5| , the self-similar character of the expan- 
sion of an initial TG gas, which significantly differs from 
the self-similar expansion of a ID TF cloud. Thus, the 
ID expansion offers a way to clearly discern between TF 
and TG regimes, and in between. We justify the valid- 
ity of the employed formalism for the expansion problem 
by comparing the hydrodynamical and the BF mapping 
results. 

We consider in the following a dilute gas of N bosons 
confined in a very elongated harmonic trap with radial 
and axial frequencies ui p and ui z (ui p 3> <jj z ). If the in- 
teraction energy per particle is smaller than the zero- 
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point energy tko p of the transversal trap, the system can 
be considered effectively as fD. We first briefly review 
the formalism introduced in Ref . |L6| . After approximat- 
ing the interparticle interaction by a delta function, the 
Hamiltonian which describes the physics of the ID gas 
becomes 
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where m is the atomic mass and (?id = — 2h 2 /mam- The 
one-dimensional scattering length is aiD = (— a p/2a)[l — 
C(a/a p )] p^] with a the three-dimensional scattering 
length, a p = yf2h/mbj p the oscillator length in the radial 
direction, and C = 1.4603 .... As shown by Lieb and Lin- 
iger (20) , Hf D can be diagonalized by using Bethe Ansatz 
p2fl . For the thermodynamic limit, a ID gas at zero tem- 
perature with a given linear density n, is characterized 
by the energy per particle 
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where 7 = 2/n|aiD|. The function e(7) fulfills 
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where g (x|7) and A (7) are the solutions of the LL system 
of equations p0| 
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We assume next that at each point z the gas is in local 
equilibrium, with local energy per particle provided by 
Eq. (^). Then, one can obtain the corresponding hydro- 
dynamic equations for the density and the atomic veloc- 
ity 
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is the Gibbs free energy per particle. Inverting the cor- 
responding Madelung transform, ip = y/nexp(iS), with 



v = (h/m)(dS/dz), one can reformulate Eqs. ([7a]) and 
( ffh| ) in the form of a nonlinear Schrodinger equation 
(NLSE) of the form: 
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Eq. (||) presents similar limitations as those of the NLSE 
of Ref. pl| . Its validity for the problem under con- 
sideration is discussed below. Note that for the case 
of n|ai£i| — > 00, one obtains 4>{n) = gmn, retrieving 
the ID Gross-Pitaevskii equation (24|, whereas for the 
case n|an)| — > 0, one gets 4>{n) = n 2 Ti 2 n 2 /2m, and 
the NLSE of Ref. j2^] is recovered. The system has 
only one control parameter [ fl6|| , namely 77 = n!p F |ai£)|, 



[(9/64)Ar 2 |aiD|/ a 4] 3 is the TF density, 



where n^p 

with a z = ^h/mui z . The regime 77 3> 1 corresponds to 
the TF limit, in which the stationary-state density pro- 
file has a parabolic form. On the other hand, the regime 
tjCl corresponds to the TG regime, which is character- 
ized by a stationary-state density profile with the form 
of a square root of a parabola. 

We have employed Eqs. (|),(|),(|),(|) and (§) to simu- 
late numerically the expansion of a ID gas when the axial 
confinement is removed, i.e. 10 z = |25| ]. In our sim- 
ulations we have employed a Crank-Nicholson method. 
Special care must be paid to the spatial and tempo- 
ral integration steps, due to the long integration times 
needed, the velocities acquired during the expansion, and 
the larger nonlinearity in comparison to the case of the 
standard GPE. 

In order to clearly understand the physics of the ex- 
pansion dynamics, let us consider the case of a bosonic 
cloud which is initially TF-like (n ~ (1 — x 2 /R 2 ), with R 
the corresponding TF radius). This would be the case of 
an initial r\ 3> 1. During the course of the expansion, the 
cloud density decreases, following in the first stages a self- 
similar TF solution n(z,t) — n(z/b(t),t — Q)/b(t) with 
b = uJ z /b 2 fj6| , p7f . However, as the density decreases, 
the gas enters from the large n|aiu| regime into the low 
n|aii}| regime. As a consequence, the functional depen- 
dence of <f>(n) changes throughout the whole cloud, and 
the expansion becomes no more TF self-similar (Fig. @). 
When this happens, the density profile departs from a 
parabolic TF profile, and asymptotically evolves towards 
a square root of a parabola shape, i.e. the density profile 
becomes TG-like. Similar behavior is observed for any 
initial value of 77. 

In our simulations we have considered due to numeri- 
cal limitations 77 values close to 1. We have analyzed in 
detail the dynamical transition from the initial density 
profile towards a TG shape, and observed that during 
the expansion the density profile presents at any time 
the form 



n(z,t) = C(t) 1 - 
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where R(t) is the radius of the cloud, and the expo- 
nent s(t) takes the value s(0) = 1 for an initial TF 
gas. The normalization constant is of the form C{i) = 
(N/^R{t))(T(s{t) + 3/2)/T(s(i) + 1)). The values of 
R(t) and s(t) are determined at any time by means of 
a nonlinear least-squares fitting algorithm. In order to 
check the validity of the fit, we have also considered the 
normalization constant as a fit parameter, and compared 
the obtained value with the expected value C(t). The 
difference is less than 0.1%. As observed in Fig. |^ (for 
rj = 1), the function s(t) decreases monotonically in time, 
and it will asymptotically reach the value 0.5 The func- 
tion s(t) presents two clear time scales. It decreases fast 
during the first stages of the expansion (few axial trap 
periods), but the final convergence towards the TG pro- 
file is significantly slower. The latter is expected since 
the exact n 2 dependence of the functional <fi(n) just ap- 
pears asymptotically. We have analyzed the expansion 
for larger initial values of 77, and observed a similar be- 
havior, although for larger values of rj the convergence 
towards s = 0.5 requires longer times (for 77 = 3 it should 
occur at around 350 axial trap periods). For practical 
purposes, if 77 is not sufficiently close to 1, only the first 
stage of the evolution of s(t) will be observable, since 
for longer time scales the density will significantly de- 
crease. For 77 ^> 1, the usual TF self-similar solution 
(s(t) = const = 1) is retrieved for any practical purposes. 
Note, however, that the transition into a TG profile dur- 
ing a ID free expansion is in principle unavoidable for 
whatever initial condition, i.e. the s coefficient tends to 
0.5. This is supported by the fact that the density profile 
fulfills at any time Eq. ( |l0| ) , and that the only self-similar 
solution supported by the NLSE (|9|) at very low densities 
is with s = 0.5. 




FIG. 1. Cloud width *J (z 2 ) as a function of time. The 
solid line is for rj = 1, (j r = 27T X 20kHz and N = 200 atoms 
(uj z = 2-7T x 1.8Hz at t = 0). The dashed (dotted) lines are 
the self-similar ID TF (TG) solutions. 

Once shown that Eq. (||) predicts that during a free ID 
expansion the density profile dynamically becomes TG- 
like, let us analyze the validity of the equation for the 
problem under consideration. As shown in Ref. (IJ], the 
hydrodynamical approach should be carefully employed, 



since it overestimates the coherence in the system. In 
order to check that Eq. (||) provides the right physical 
picture in our problem, we have calculated the free ex- 
pansion of an initial TG gas using both the BF map, and 
the NLSE. 

From the BF map one obtains that the dynamics of 
the density profile for an impenetrable gas of bosons is 
given by p) 
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where (j) n (z,t) denotes the time-dependent wavefunction 
of the n-th eigenmode of the original axial harmonic os- 
cillator. The expansion dynamics for each <f> n is obtained 
analytically by means of the corresponding Green func- 
tion in free space 

G(X - Z', t) = -i( V2 e im{z-z'f/2Ht_ (12) 

From Eq. (|ll]) one obtains a self-similar solution of the 
form 
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Note, that for times t 3> l/ui z the scaling coefficient 
i/l + iv 2 t 2 becomes cu z t, whereas for the case of a ID 
TF self-similar solution the scaling coefficient becomes 
\/2uj z t. Consequently, the expansion of an initial TG 
and TF gas is significantly different. This property could 
be employed to discern between the two regimes. 



time (l/(O z ) 



FIG. 2. Time evolution of the exponent s(t) for 77 = 1, for 
the same parameters as in Fig. |l|. The exponent decreases 
monotonically, and will asymptotically approach 0.5. 



From the corresponding hydrodynamic equations QTe 
and (ffb]), one can easily prove that the same self-similar 
solution dl3| ) for the density is obtained from Eq. (|]) 
in the limit of n|ai£)| — * 0, i.e. using the equation of 
Ref. pi] . Therefore, Eq. (|J) accurately describes the 
expansion dynamics even for the extreme case of a TG 
gas, and it is thus expected to describe well the expansion 
for intermediate regimes between the TG and the TF 
limits, where the coherence is not yet completely lost. 
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In this Letter, we have studied the dynamical tran- 
sition from a quasi-lD BEC into a TG gas during the 
expansion. Our analysis is based on a NLSE with vari- 
able nonlinearity, which generalizes for arbitrary interac- 
tion the extremal cases provided by the Gross-Pitaevskii 
equation (n|aic| — * oo) and the equation of Ref. 
(n|an)| — * 0). We have shown that even if the initial 
cloud possesses a TF profile, the density profile acquires 
asymptotically a TG shape. We have analyzed in de- 
tail this transition, and characterized the shape of the 
cloud in the intermediate stages. We have evaluated by 
means of a BF map the exact expansion dynamics of 
a TG gas and shown that the expansion is self-similar 
with a significantly different scaling law compared to a 
TF gas. We have additionally shown that the NLSE ap- 
proach provides exactly the same self-similar solution as 
the BF map for the case of a TG gas, and it is therefore 
expected to describe well the expansion for any interme- 
diate regime. 

Let us additionally point out that the NLSE (||) also 
provides the excitation spectrum of the ID Bose gas in 
intermediate regimes between TF and TG, by considering 
a small perturbation around the ground state solution 
-00 (z) of Eq. (j^ ip(z) = ip (z) + Sil;(z), where Sip is given 



by 8ip{z) = u(z)e- luJt + v( Z y 



Inserting this Ansatz 



into Eq. (j|) leads to the corresponding Bogoliubov-de 
Gennes equations 



Cu(z) + no<fi' (no)v(z) = Tiuju(z) 
-Lv{z) — riQip' (jiq)u(z) = fru>v(z) 
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where 
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d(f>/dn, and C 



-(hy2m)(d 2 /dz 2 ) + mu> 2 z 2 /2 + cf>(n ) + n o 0'(n o ) - fi, 
with /i the chemical potential fixed by the normalization 
of no. Eqs. ( |l4| ) and ( |l5| ) describe the crossover from the 
TF to the TG regime for all excitation frequencies. In 
particular, we have obtained that these equations provide 
the same results as in Ref. J2^] for the lowest compres- 
sional mode. 

Summarizing, the ID expansion dynamics constitutes 
an experimentally accessible tool to discern between the 
different interaction regimes in a ID gas, and addition- 
ally provides a way to dynamically accomplish the TG 
gas. Unfortunately, the method employed in this Let- 
ter does not allow to analyze the fundamental problem 
of decoherence when entering the TG regime. The solu- 
tion of this problem requires to extend the exact results 
of Rcfs. 0Jl| to the case of inhomogeneous time de- 
pendent Bose gases with finite interactions, in which the 
BF mapping is not exact. Such analysis is beyond the 
scope of this Letter and it will be the subject of future 
investigations. 
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